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We derive a formula describing the evolution of tip-splittings of Saffman- Taylor fingers in a 
Hele-Shaw cell, at zero surface tension. 
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Fingered patterns characterize the non-equilibrium growth processes of many systems, 
e.g. the dendritic shape of snowflakes, the outline of bacteria colonies growing in stressed en- 
vironments, electro-deposition, dielectric breakdown, and viscous fingering formed by forcing 
a non-viscous fluid into the center of a Hele-Shaw cell filled with a viscous fluid (Saffman- 
Taylor flow)|lJ. Repeated events of tip-splittings and side branching during the growing 
process, manifest themselves in complicated fractal patterns. Understanding the structure 
of such patterns is a challenge for theorists. 

When the growing patterns are self similar, one expects that their global structure can 
be deduced from the basic elements of the growing process, such as tip-splitting and side- 
branching. Indeed, studies of theoretical models of diffusion limited aggregation^, and their 
generalizations (e.g. the dielectric breakdown model) demonstrate that the fractal dimension 
of the corresponding patterns is related to characteristics of the tip-splitting events 0] (side- 
branching in these models is negligible). 

In this work we focus our attention on the fingered patterns generated by Saffman- Taylor 
(ST) flows in a Hele-Shaw cell, in the radial geometry. In particular, we shall be interested in 
the evolution of tip-splitting in the singular limit where the surface tension associated with 
the interface between the viscous and the non-viscous fluids approaches zero (or alternatively 
high pumping rate of the non- viscous fluid). The ST dynamics in this limit, which we shall 
refer to as the idealized ST problem, has been shown to be integrable[4]. Therefore, it 
constitutes an important paradigm in the field of non-equilibrium growth processes. 

The idealized ST problem admits a large class of exact solutions j^, f], Q| • Among t 
one may find solutions which resemble various forms of tips-split at large surface tension ^ 
These solutions are characterized by a smooth evolution of the bubble counter for any time. 
Yet, these solutions lack universality since they depend on precise details of the initial 
conditions. 

Nevertheless, there is another class of solutions of the idealized ST problem: Solution 
which exhibit a cusp-like singularity after finite time. A generic ST bubble always develops 
a cusp. The cusp singularity dominates the evolution in its vicinity, and therefore it is 
expected to be universal. Thus, our aim is to describe the local dynamics near cusps of 
tip-splitting. 

Our central result is, thus, a formula which describes the evolution of tip splitting in 
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FIG. 1: The evolution of tip splittings in the Saffman- Taylor problem at zero surface tension. The 
contours represent snapshots of the evolution as function of the time. The evolution given in (^Q), 
is characterized by one parameter, eft, which controls the asymmetry between the two generated 
fingers. Panel A shows the symmetric evolution, while panel B represents a typical asymmetric 
behavior. 

time. It has the form 

z( S ,t) = s 2 +u,(t)+ v ;l t] . a) 

W^t) + IS 

where z = x + iy is a complex coordinate on the bubble contour, — oo < s < oo parametrize 
the curve, and u^it), v^(t) and w<f,(t) are functions of the time, t. These functions which will 
be calculated in what follows, depend on a single parameter, 0, governing the asymmetric 
shape of the evolution. In Fig. 1 we depict contours obtained from ((TJ which represent 
snapshots of the tip-splitting evolution as function of the time. 

A problem one encounters when trying to describe the tip-splitting evolution is that a 
naive extension of the idealized ST dynamics beyond the cusp singularity is impossible. 
One may try to compute the ST dynamics, in the vicinity of a cusp, by introducing some 
infinitesimal surface tension. This approach indeed regularizes the problem but the resulting 
equations are nonintegrable, and the analytical treatment becomes very complicated |9|. 
There are some exact steady-state solutions, with surface tension, but these do not exhibit 
tip-splitting jlpl |. 

Alternativel y, t he cusp singularity may be resolved by a method known as "dispersive 
regularization"|ri| [l^. which preserves the integrable structure of the problem. This ap- 
proach is similar, in spirit, to the construction of Gurevitch-Pitaevskii solutions for the KdV 
equation Its application to the ST problem is discussed extensively in 17]. Here we 



shall focus on time intervals where these aspects of the regularization are irrelevant. 



We begin by recalling the ST problem. The local velocity of a viscous fluid in a thin cell 
is proportional to the gradient of the pressure v oc Vp, where p(z) is a harmonic function 
due to incompressibility. In the zero surface tension limit, the pressure may be taken to be 
equal to zero on the perimeter of the bubble, while at infinity it diverges logarithmically. At 
constant flow rate, the area of the bubble grows linearly with time, t. We set the flow rate 
such that the area is nt. The other parameters, viscosity and width of the cell, are chosen 
such that the pressure is equal to the velocity potential. 

An important feature of the ST dynamics, at zero surface tension, is the conservation of 
the harmonic moments: -4^ = 0, where 
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tk = ~ I . f?z~ k , fe = l,2,.... (2) 
nk jvisc. fluid 

and the integration is over the domain occupied by the viscous fluid. This property implies 
that the idealized ST dynamics is integrable, i.e. the bubble's contour can be determined 
from the set of harmonic moments, t^'s, and the area nt. Our goal is to describe the evolution 
of this contour in the vicinity of the tip-splitting. 

In order to do, so it is instructive to introduce a potential defined as V(z) = t\og(z) + 
YjktkZ k ■ I n particular, it will be convenient to work with sets of harmonic moments for 
which the potential can be resummed as 

N 

V(z)=t\og(z)+t 1 z + J2^og(z-q l ). (3) 

i=l 

Thus, the set {tk}™ =1 is replaced by a new set of parameters {fa}^ =1 and {<&} i=1 known as 
Miwa variables in the soliton literature^]. In particular and fa are, respectively, the 
location and the weight (assumed to be real) of the z-th Miwa variable. 

The potential defined above encodes the shape of the bubble in a rather complicated 
way. A function from which it is easier to extract this information is the Schwarz function. 
This function has proved to be a useful tool for the analysis of the ST problem (see e.g. 
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10[[15[|16||)- A Schwarz function, S[z), of a given contour C satisfies the relation S(z) = z 
(where z is the complex conjugate of z) for points z which lie on the contour C. Away 
from the contour it is defined via analytic continuation. The simplest example is a circular 
contour of area 7rt, whose Schwarz function can be easily deduced to be S(z) = -. 

The Schwarz function is related to the set of harmonic moment by the contour integral 
tfc = §c S(z)z k dz. It implies that, on the exterior of the bubble, S(z) has the same singular 



structure as (as can be seen by deforming the contour integral such that it envelopes 
the singularities of S(z) on the exterior of the bubble). Thus the analytic structure of S(z) 
in the exterior domain, can be extracted from the potential, V(z). 

To reveal the analytic behavior of S(z) in the interior of the bubble one may employ an 
important property of the Schwarz function known as the unitarity condition, 

S(S(z)) = z, (4) 

where the complex conjugation of a function, f(z), is defined as f(z) = f(z). The equation 
above trivially holds for z G C and by analytic continuation all over the complex plane. 
S(z) maps points from the exterior of the bubble (where the analytic structure is known) 
to points on the interior domain. Thus all singularities of S(z) can be identified, and hence 
its precise analytic form, as function of the Miwa variables and the time. 

We proceed by assuming that the Schwarz function is an algebraic function, i.e. a function 
defined on some algebraic Riemann surface. Let us explain what is an algebraic Riemann 
surface, and then describe how to define a Schwarz function on this surface. Consider a 
polynomial equation of the form P(z, f(z)) = 0, where P is a polynomial in two variables. 
Solving this equation we obtain for each z a set of solutions f(z). In general f(z) will have 
branch cuts. A simple example of such a polynomial is P(z, f) = f 2 — Hf =1 (z — Aj), where 
we obtain f(z) = \Jllf =1 (z — A«). In general there will be n solutions of / for each z. Thus 
one may introduce n copies of the complex plane, where each copy is associated with a well 
defined value of f(z). Clearly, on each copy of the complex plane f(z) is discontinuous along 
the branch cuts. To define it as a continuous function we may paste together the various 
sheets of the complex plane along the branch cuts. Then, together with a choice of local 
coordinate systems around each point one obtains a Riemann surface composed of n sheets. 

Since we assume that the Schwarz function may be defined on such a Riemann surface, 
one should assign for each point on the Riemann surface, a value of Schwarz function, S. 
In order to specify such a point we must indicate the sheet-index i (where 1 < % < n), and 
a complex number z. The index, i, will specify on which copy of the complex plane the 
point lies, and a complex number, z, will specify the coordinate on that copy. There will be 
one copy of the complex plane, which will be termed as "the physical sheet" , on which the 
bubble lies. On this Riemann sheet S(z) will be equal to z on the perimeter of the bubble. 
All other "unphysical" Riemann sheets will be considered as parts of the "interior domain" 
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since there are no branch cuts on the exterior domain (where the pressure is harmonic). 

Let us now consider a situation where the potential has the form ©, and show that the 
number of unphysical sheets is the number of Miwa variables. On the exterior of the bubble 
S(z) has the following singular structure determined by 



For the point q iy unitarity (jJJ) implies that S(S(qi)) = S(oo) = q; t . Thus a point at infinity 
should be mapped to the point cji. This point should be located on one of the unphysical 
sheets since on the physical sheet S(oo) = t±. Conversely, the value of the Schwarz function 
at infinity on an unphysical sheet correspond to a particular Miwa variable. From here on 
we consider the case of two Miwa variables, which is the minimal model exhibiting a tip 
splitting scenario. The corresponding Riemann surface is thus composed of three sheets. 

For actual calculation of the contour's evolution it will be convenient to use a conformal 
mapping which maps the exterior of some "source domain" (in £-plane) to the exterior of 
the physical bubble, the "target domain" (in z-space). The source domain is usually taken 
to be the unit circle, however, we found it more convenient to use a bubble with a cusp. The 
advantage in using this source domain is that we can choose the mapping to be non-trivial 
only around the tip, while everywhere else it would be approximately proportional to the 
identity map. 

We will study the evolution near the tip in the case where the weight of one of the Miwa 
variables, is small compared to the other. This suggests choosing the source domain to be 
the one Miwa-variable bubble at the point where a cusp is formed. The Riemann surface 
in this case is composed of two sheets as described above. If we assume this surface to be 
of genus zero, then by the Riemann-Hurwitz theorem, the number of branch points is two. 
Thus the Riemann surface is given by the polynomial equation R 2 = (£ — Ai)(£ — A2). A 
function defined on this surface (the Schwarz function in particular) contains a branch cut 
which extends from point Ai to the point A2. An example of such a function is R(C)- 

The 1-Miwa bubble is characterized by four parameters ti,t,q and /1, since the location 
of the branch points depends on these parameters, one can calculate the area of the bubble 
and its first harmonic moment, £1, from the location of the branch points Ai and A2. We 
will fix fi = —q = 1 and treat Ai and A2 as the parameters which describe the bubble. If we 
also demand that the bubble is at the moment where the cusp forms, we can characterize 
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the bubble by a single parameter, say Ai. We assume that Ai < A2 and that Ai is of order 
5, where S <C 1. This assumption implies that the global bubble shape is dominated by the 
cusp. Then the bubble area is of order 5 3 while the first harmonic moment, t\, is 1 + 0(5 2 ). 
The Schwarz function, <j((), of such a bubble is: 

2«(C) = ( il -l) + F| - — + ( ( - tl )(C + l) ' 

as can be ascertained by examining this function's singular structure and that it satisfies 
the unitarity condition (j3J). The fact that this function describes a bubble with a cusp can 
be checked by considering the behavior of the solution of a(() = ( near Ai (where the cusp 
is located). 

Up till now we have characterized the "source domain". We would like, now, to specify 
the physical bubble, or the "target domain" associated with two Miwa variables. For this 
purpose we must give the mapping between the source and target domains. This mapping 
is taken to be: 

where ci, c 2 , a, (3, and 9 are parameters of the mapping and R(() = VC — ^iVC — ^2- This 
mapping can be considered as a mapping from the 1-Miwa Riemann surface to the 2-Miwa 
Riemann surface. The mapping has singularities at the infinities on each of the sheets of 
the Riemann surface associated with the source domain, which are mapped to infinites on 
different sheets on the target domain, as well as at the point 9 on the unphysical sheet. This 
point, 9, is mapped to an infinity on a third sheet of the target Riemann surface. Thus (0) 
indeed maps a two-sheeted Riemann surface, which is associated with a 1-Miwa bubble, to 
a three sheeted Riemann surface. 

That, indeed, the target Riemann surface is associated with a 2-Miwa variable bubble, 
can be deduced from the singular structure of the Schwarz function of the target domain. 
The latter is given by S(z) = z (a (C (z))), where ((z) is the inverse map of z(C)> and <j(() 
is the Schwarz function of the source domain. 

Having S(z) one may extract all constants of motion (Miwa variables, qi, and q<i ; Miwa 
weights, /ii and ^2] and £1), and the area t, from its singular structure. These will be 
expressed as functions of the parameters of the mapping © and Ai. Solving these relations 
one may express the parameters of the mapping as functions of the time and thus obtain the 
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evolution of the contour. We take fi\ = —q\ = 1 by fixing c\ and 02- Thus the parameters 
whose time evolution is to be determined are a, (3, 6 and Ai. 

To obtain reasonably simple equations, we expand all quantities in orders of 5 and take 
the leading order. Let us assume that for some initial moment, around the formation 
of the tip, Ai assumes the value A, to leading order in 5. We now make the following scaling 
ansatz: a ~ 5 4 , (3 ~ 5 3 , v = Ai — 9 ~ 5 3 and SX = Ai — A ~ 5 3 . Then the equations we 
obtain, to leading order in 5, for the constants of motion q2 — A and fj,2 are: 



(6) 



^-2^^[l- — ) (7) 




Let us now define T\ — t\ — (l + 3A 2 — 7A 3 + yA 4 ^ (the difference between the first 
harmonic moment of the target bubble and a source bubble with Ai = A to order 0(S 5 )), and 
similarly T = t — (— 4A 3 + 18A 4 — 63A 5 ). With these definitions, analysis of the singularities 
of S(z) yields: 

T = 2^{a)X + A(3X 2 -12X 2 5X, (8) 
T x = R(a) - 2/3A + 6A.5A. (9) 

The solution of Eqs. (013) give the parameters a, (3, v and 5X in terms of q2, fJ-2, T and Ti, 
which define, in turn, the conformal mapping of the contour on the target space as function 
of the time. To write down the solution of these equations it will be convenient to define a 
shifted rescaled time 5t = — ^j^p = and introduce two functions £ and rj which satisfy the 
nonlinear equations: 



(y/25t 5t 2 \ . ^ 1 
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where (j) = S(-2=)) is the asymmetry parameter. The solution of these equations gives rj 
and £ as functions of the rescaled time, 8t, and (p. With the help of these functions we may 
write the solution of Eqs. (JHIEJ) as: 



-X V R 



v = — g— +«V^' 
4« = 3 Sfe )-3A-5(,-fl-^| + g-| 
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Given this time dependence of the parameters of the conformal mapping we are in a 
position to describe the contour dynamics in the vicinity of the tip-splitting. For this purpose 
it is sufficient to focus on the image of the source domain around the cusp. The shape of the 
source domain near the cusp is given by the universal form y = ASx 3 ^ 2 , were 5x = x — Xi, and 
A is some constant. Close enough to the cusp y <C 8x, and therefore we may assume that 
y ~ (this assumption can be proved to be consistent with the expansion in the parameter 
5 performed above). Thus one has to find the image of the ray x > Ai under the mapping 
z(() to leading order in 5. The result is given by Eq. (JTJ), where the functions u^t), v^t), 
and Wtpit) are: 

MO = ~ o ~ " 



2£' 

v+(t) =6t(l- , 

The above equations together with (0) describe the evolution of a tip-splitting of the ST 
problem at zero surface tension (Fig. 1). The form of the tip-splitting depends on a single 
parameter, <f>, which controls the asymmetric shape of the evolution. Since our derivation 
of the tip-splitting formula makes use only of local properties near the tip, it is suggestive 
that this evolution, for short times, is universal. Namely, the tip splitting evolution is 
characterized by one parameter, cf>, independent of the shape of the bubble on large scales. 
The shortcomings of our analysis is that it does not include the influence of surface tension. 
Therefore it breaks down after a short time due to the formation of cusps. The method 
of dispersive regularization can be employed again to continue the evolution. Thus, in 
comparing our theoretical prediction with experiment, one can hope for agreement only 
within a limited time interval near the initial stage of the tip-splitting. 
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